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^ i ABSTRACT 

We study a theory which generalizes the nonminimal coupling of matter to gravity 
by including derivative couplings. This leads to several interesting new dynam- 
ical phenomena in cosmology. In particular, the range of parameters in which 
inflationary attractors exist is greatly expanded. We also numerically integrate 
the field equations and draw the phase space of the model in second order approx- 
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O . imation. The model introduced here may display different inflationary epochs, 

generating a non-scale-invariant fluctuation spectrum without the need of two 
^ ■ or more fields. Finally, we comment on the bubble spectrum arising during a 

Q\ ■ first-order phase transition occurring in our model. 

"o 
cr 

u 

1 Introduction 

^ . Scalar fields in General Relativity has been a topic of great interest in the latest years, 
mainly because the scalar field dynamics allows to investigate the detailed features of the 
early Universe. Without the need of a specific equation of state, a Universe filled by a 
phenomenological scalar field leads to an accelerated phase of expansion, the inflation, in 
almost any kind of self-interaction potential. As it has been shown in chaotic models of 
inflation0, the accelerated phase, either power- law! or quasi-exponential, is a phase-space 
attractor for most of the initial conditions, and this result also extends to a large class 
of inhomogeneous and anisotropic space-times!. The class of successful models has been 
enlarged to scalar fields with a nonminimal coupling (NMC) to the curvature scalar R in 
the gravity Lagrangianl, commonly in the form of a term y/—gf{(j))R, where /(</>) is a 
function of the scalar field <fi. The motivations for this step are manyfold: the idea that the 
fundamental constants are not constant, the "machian" theory of gravitation embodied in 
the Jordan-Brans-Dicke theory!, the renormalizing term arising in quantum field theory in 
curved spaced, the possibility to have gravity as a spontaneous symmetry-breaking effect!, 
the Kaluza-Klein compactification schemei, the low-energy limit of the superstring theory!. 
Moreover, the NMC term has been employed to produce an oscillating Universe,0 to reconcile 
cosmic strings production with inflation Jill to generate a modified Newtonian dynamics able 
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to model flat rotation curves in galaxies.El Last but not least, an NMC term allows to solve 
the "graceful exit"_problem of the old inflation by slowing the false-vacuum expansion, as in 
extended inflation^. Recently, the NMC models have been generalized in various directions: 
models with different coupling functions^, with generalized coupling to the inflaton sectoral, 
with fourth-order gravity^, and with a coupling to a dark matter sectored. 

In this Letter we wish to further explore the influence of nonminimal couplings in cosmol- 
ogy by introducing the nonminimal derivative coupling (NMDC) to gravity. In this class of 
models, the coupling function is also function of the derivatives of 0: / = f((f>, <f>. tli , 4> ;flu , ...). 
Derivative couplings, although rares, are not a novelty in physics: the field theory of scalar 
quantum electro-dynamics includes derivative couplings between the electro-magnetic vector 

and the scalar field 0; this kind of interaction is indeed required by the U(l) gauge- 
invariance of the theory. In addition, if the idea underlying the nonminimal coupling theory 
is that the Newton constant G itself depends on the gravitational field source mass, it seems 
more natural to couple the curvature to the energy-momentum tensor of the matter, intro- 
ducing the terms TR and T^ U R' 1U , which in fact contain derivative couplings. The NMDC 
model can be thought of as the scalar field formulation of the hypothesis that G = G(p), 
where p is the energy density of the gravitational field source. 

One of the aims of this Letter is to see whether the presence of NMDCs allows to find 
inflationary attractors in models where otherwise they are not present. We studied a large 
class of NMC models in Ref. ([0|). For future reference, let us summarize the main results 
of Ref. (|T8|). Hereinafter, we will use / = f(4>) referring to the non- derivative terms coupled 
to the curvature scalar, and f\ = fi{<j),<j) ;a ) when we refer to all the terms, including the 
derivative terms, coupled to the curvature scalar. In Ref. (fT3) the coupling function f(<f>) 
is left as general as possible, and the potential is assumed to be V = Xf M - It is also 
assumed that /(</>) is semi-positive definite and that, for large \<fi\, it grows monotonically 
faster than <fi 2 . The case / ~ <ft 2 is the standard choice, and has been extensively treated 
in the literaturei. The relation between the potential and the coupling function assumed 
in Ref. (|18D allows to determine the class of models with successful inflationary attractors; 
here "successful" means that the inflation eventually stops and a Friedmann stage begins. 
In Ref. ([I%D and in the present work the inflationary attractors are all asymptotic, in the 
limit of large \<f>\; in this limit the influence of the coupling is clearly larger. Moreover, 
the spirit of chaotic inflation is that whatever the initial conditions of the Universe are, the 
dynamical trajectory falls onto an inflationary attractor. This is verified only if the attractor 
itself extends to very large values of the field variables, until the classical description loses 
meaning (i.e. near the Planck boundary). In Ref. (|18|) it is shown that such asymptotic 
attractors exist only if 2 < M < 2 + y/3, whatever functional form of /(</>) is considered. 
In this range, it is found for the cosmic scale factor a(t) a power-law inflation a ~ t p with 
exponent 

P = 3 + (M-2)3 

P (M-2) 2 (M- 1) 1 J 

In the conformally rescaled frame the power-law exponent has a much simpler form, p = 
3/(M — 2) 2 . If M = 2 a deSitter exponential inflation is found. Further, in the cases in 
which a successful inflation exists, the spectrum of primordial fluctuations is independent of 
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/(</>) and depends on the wavenumber k according to the power- law fcV(i-p). Finally, in Ref. 
(pp it is shown that the fluctuations have a Gaussian distribution. In the present work we 
extend some of the cited results to derivative couplings. We assume as in Ref. (|18|) that 
V = \f M , but now we will freeze the functional form of /(</>) to a power-law, /(</>) = 4> 2m . 
We will show that the narrow range in which the parameter M is confined in NMC models 
is greatly expanded when one includes a derivative coupling. 

It rises natural the question of whether the conformal transformation!!!, which allows one 
to put complicated unconventional gravity theory in pure Einsteinian form with one or more 
scalar fields, works also when derivative couplings are present. We will show schematically 
that a conformal rescaling g^ u = e 2 ^ g^v cannot recast our theory in Einsteinian form. In order 
to recover the Einstein field equations the metric transformation should be generalized to a 
Legendre transformation, as in Ref. fl2"0|). This would introduce in the equations additional 
tensor fields, instead of additional scalar fields as in non-derivative cases. Due to this reason, 
we avoided here the use of metric transformations. 

The plan of this paper is as follows: in the next section we introduce the derivative 
couplings, in Sect. 3 we derive the equations of motion for the fields, in Sect. 4 we perform 
the second-order approximation, in Sect. 5 we discuss the phase space of some models and 
finally in Sect. 6 we present the conclusions. 



Derivative couplings 



Terms with NMDCs generalize the Einstein-Hilbert Lagrangian. The most general gravity 
Lagrangian linear in the curvature scalar R, quadratic in <ft, containing terms with four 
derivatives includes all of the following terms: 



■a/3 ) 



L 5 = vcj)(j). a R' a ; Lq = a< 



9UR, 



(2) 



The constants /z, r, ... have the same dimensions in natural units as the Newton constant G, 

namely mass" 2 . Here and in the following the conventions are: signature (H ) ; units 

87rG = c = 1 . Due to the following relations 
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three of the six terms L\ — Lq can be neglected (also taking into account the Bianchi iden- 
tities). Other total divergences, like (cf) 2 R^) tll and (0 2 i?^) ;At!/ are linear combinations of 
Eqs. (^). Nevertheless, the equations remain hopelessly complicated. We take here a first 
step. We keep only the derivative term L±, together with the NMC correction Lq = £f(<j))R. 
This particular choice turns out to be the simplest one. Moreover, the terms studied here 
are those appearing in the simplest "natural" coupling, namely TR. In this way we may 
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draw an useful comparison between Ref. (|I5D , where Lq was the only correction to Einstein 
gravity, and the present work. We are then left with the following Lagrangian density 

fl&ha) = (4) 

3 Field equations 

The set of equations derived from (H) is 



G>„(1 - fx) = Q^fl - fl^u + H<P;v<P;vR + T t 



= </> itl <f> iV - ^g^<P; a (p ]a + V((j))g^ u , (5) 

while the scalar field equation is 

□0(1 + fiR) + fi(f); a R a - if'R/2 + V = . (6) 

It can already be seen that, even if the new terms introduced by the derivative coupling (let 
us call them //-terms) do not determine the critical points of the theory, they do contribute 
to the asymptotic dynamical properties of Eq. (^|). Indeed we will see that in the limit of 
large R and R. a (in other words, in the early Universe), the //-terms will sensibly modify the 
behavior of the non-derivative theory. 
The trace of Eqs. (|5]) reads simply 

fl(^/-l)=T + 3D/ 1 , (7) 

where T = —(p- a (p' ,a + 4V. The derivative of the Ricci scalar, to be used in Eq. @, is 

In a Friedmann- Robert son- Walker (FRW) spatially flat metric with scale factor a(t), the 
(0, 0) component of ([5]) is 

H\l - fx) = Hh + + Uh* + V) , (9) 



3 V2 

where H = a/a and where our conventions are such that R = —6H — 12H 2 . From now on, 
we confine ourselves to a FRW spatially flat metric. The system (||-|9D is closed. We have 
five degrees of freedom H, <p, <ft, <f), <ft in and indeed we have one equations of fourth order in 
<p [Eq. ([])] and one of first order in H [Eq. (J|)]. The solution of this system appears very 
difficult. In the next section we will make use of a second-order approximation to derive the 
asymptotic properties of our model and to sketch its phase space. 

Let us conclude this section showing how a conformal transformation operates on the 
NMDC theory. Ordinary NMC theories are recast in Einstein form by the metric rescaling0 
g^u = e^g^, with 

e 2w = \dL/8R\ . (10) 
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Adopting the same formula, in our case we get e 2uj = 
tensor G^ u = R^ u — Rg„ u /2 trasforms according to 



IWi 



Now, we have that the Einstein 

(11) 



Inserting (|i~I|) into (||), with the choice CiTOf) , one has 
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3e 



hwhv - ^9^h;J[ + lie 



-2u . 



'\liV\v 



R + e~ 2 "T t 



/J.U 1 



(12) 



which is not in Einsteinian form. Any other choice of the conformal factor would not cancel 
the second derivatives of f\ which are present in Then, the right-hand-side of ( |T2"D could 
not be written as a scalar field energy-momentum tensor. Notice that R at the right-hand- 
side is not trasformed. If one also trasforms R in (0), additional terms like (Da;)0 ;A 
arise. 



4 Second-order approximation 



In slow-rolling motion one neglects all terms with more than one time derivative in the 
field equations. In this section we are more general, and we examine the field equations by 
neglecting all terms of order higher than the second one. This gives us three advantages: the 
equations are notably simplified, we can calculate analytically some inflationary solutions, 
and we can plot the phase-space portrait of the model. Neglecting all terms higher than 
second order in (B) we have 



□0(1 + nRW) + - if'R {2) /2 + V = 



(13) 



The indexes in parentheses denote the order to which R and R are to be calculated. From 
Eqs. (|7|,!D an d for £/ ^> 1 we have 
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We now explicitely adopt the functional relation V 
Eq. fljjD then writes 



(14) 
(15) 

(16) 

Xf M . This will simplify our equations. 
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f M ~\M 



(17) 



Let us call A the coefficient of □</> and B the coefficient of f'(f) 2 . In order to avoid singularities 
in Eq. (\L7\) one should have a nonvanishing coefficient A for £/ 3> 1 (in the opposite limit 
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the equations reduce to the minimally coupled model and the regularity follows automat- 
ically). A sufficient condition, although not a necessary one, is that £, /x < 0, along with 
the already imposed condition / > 0. This assumption greatly simplifies the investigation 
of the dynamical properties of the model under study and will be adopted in the following. 
It also makes clearer the comparison with previous work in ordinary NMC theories, where 
the condition £ < is often adopted to get a regular metric rescaling. Notice that £, \i < 
ensures the positivity of the effective Newton constant in the Lagrangian density ([§). 

The d'Alambertian operator reads □</> = + 3if in FRW metric, where the Hubble 
function H is to be taken from Eq. (H) (to the first order): 

- 2 , 1 fl l2 



The asymptotic behavior of Eq. (|T7|) can be now easily discussed in the limit in which 
some of the terms in A and B can be neglected. A first consideration is however immediate: 
for M < 2 the second derivative of the effective potential in Eq. (|I7]) is negative, in the 
hypothesis that / grows monotonically with \(f)\. As a consequence, any attractor solution of 
Eq. (|H]) leads to ever growing /; the Friedmann behavior cannot be reached. Let us make 
an example. If we assume f(<p) — > <p 2m for large \(f)\, it turns out that if M < 2 — 1/m one 
has that f' 2 dominates in A and /" in B. Then Eq. (|D]) simplifies to 

□ / _[2A(M-2)/3£]/ M = (19) 

(where □/ = f'n</> + f"<j) 2 ). In this case the influence of the /i-terms disappears; the model 
reduces to the ordinary NMC theory, already discussed in Ref. ([18]). There it was found that 
for M < 2 no successful inflation was allowed. Indeed, in Eq. ( |19"1) the effective potential 
derivative is negative definite (since £ < and M < 2). As already discussed, the attractor 
solutions of Eq. fll~9|) are directed toward increasing /; the Friedmann behavior will never be 
recovered. The same conclusion holds for any f(4>), provided that f' 2 ^> f M , f. 

In the special case M = 2 one sees that <p = <ft = is a solution of ([H])- One has then 
an asymptotic deSitter inflation, analogously to what occurs in NMC models. We will not 
give details of this case, since the dynamics depends on the behavior of V and / at small <ft, 
while here we are interested mainly in the asymptotic properties. 

Let us come to the third, more interesting, case, namely M > 2. Now f M dominates in 
A and / M_1 in B. Eq. (|17|) then reads 

□0 + 40 2 (M-l) + -^f(M-2) = O. (20) 

It worths remarking that the effective potential is U(4>) = £(M — 2)(4/i) -1 /(</>); the coupling 
function f(<fi) reveals itself as an effective potential for the theory. Notice also that, contrary 
to the previous case, the constant /z plays here an important role. Let us put now /(0) = x m , 
where x = cf 2 . It then follows 

x 2 £N 

Ux [1-2N -2m] + —x m = 0, (21) 

2x ji 
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where N = m(M — 2), for x > 0. To have a positive-definite potential derivative is sufficient 
to have M > 2 (provided that £/n>0) . 

To find explicit solutions to Eq. ( pip we need also H = H(cj),<p). From Eq. (|l"8|). 
neglecting all terms higher than first order and in the limit £/ ^> 1 , we get for M > 2 — 1 jm 
(i.e. AT > -1) 



[H 



(1)12 



A 



JV-l 



^x 2 + x m+1 



(22) 



Let us put a trial solution x = bx p in (pi]) and (|22f) . In the limit of large x we find the 
attractor 



P 
6 



iV)/2, 



(23) 



It can be easily proved that this really is an attractor; a graphic evidence is provided by 
our phase-space portraits, Figs (|T],0). Along the solution ( P3"D the term 'x is negligible; this 
justifies the second order approximation previously performed. Integrating x = bx p one gets 



x = x (l + t/r) 1 ^ 1 - p) 



(24) 



where r _1 = 6(1 — p)xq~ and x(t = 0) = Xq. For p < 1 the time constant r is negative; 
clearly our approximations break down for t — *■ r. The inflationary regime we will find also 
breaks for t — > r. Along the attractor we have the Hubble function 



h = {\/m) ll2 x {N+m)n 



and the cosmic scale factor 



a{t) = a exp [e [l - (1 + t/r) k ] } , 



where 



E 



3(|/V£|K 
N(N + 1) 



JV+l 



k 



2(N+1) 
2 + N -m 



(25) 



(26) 



(27) 



The scale factor expansion is always inflationary for t < r. For t<r the expansion follows 
the deSitter law a ~ exp[(E\k/r\)t], (k/r is always negative). The constant E sets the total 
number of e-foldings of the inflationary stage. If, for instance, £ and \x are approximatively 



1/2 

equal, values of (/>o = x around unity in Planck units are required to have E > 60, the same 
requirement commonly found in chaotic inflation. This kind of behavior is often called quasi- 
deSitter inflation. The narrow range in which viable inflation is found in NMC theories^, 
2 < M < 2 + v3, is here expanded to M > 2, without upper bound. The effect of including 
a derivative coupling is then to substantially enlarge the class of viable inflationary models. 
This is the main result of this paper. 
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It can be checked that H(t) and x(t) decrease as the time increases. As \<p\ decreases, the 
solution crosses to the NMC region, where the derivative coupling is negligible compared to 
the non-derivative one. In particular, the attractor enters the NMC region for 

* < * = [hxl ■ (28) 

As |0| further decreases, becomes smaller than unity for x < x 2 = £ _1 / m and the 

attractor enters the minimally coupled central region. For instance, in a model in which 
| /x | A = 3/2 and M — 3,m — 1, one has X\ — £ and x 2 = l/£. If x 2 > X\ the NMC phase 
disappears. 

In Figs. (0,0) we show the numerical phase space of the model for some parameter values. 
The plots are obtained integrating Eq. (fP7D for several initial conditions, and then performing 
a Poincare projection onto the unitary circle. For any trajectory we can identify four stages, 
or cosmological epochs. All trajectories start ideally from the initial singularity at <p, <fi — ► oo, 
and they fall after a short transient on the inflationary attractors. In this first stage the 
higher-order terms that we are neglecting in this section can be important. The second 
stage is represented by the outer part of the attractors, where the /x-terms are dominant; 
now the trajectories follow Eq. (|2"3"D, and the cosmic expansion is given by (EBJ). The third 
stage begins when x < x±, where the attractors enter the NMC region. Now the behavior 
will be as described in the introduction, with the power-law expansion ([!]) (inflationary if 
2 < M < 2 + v3)- Finally, the coupling terms become negligible (for x < x%) and the theory 
reduces to ordinary gravity. The details of this last stage depend on the specific form of 
V((f)) for small <ft, but eventually (perhaps after a last quasi-deSitter inflationary episode) 
the mass term will dominate, and a series of damped oscillations around the potential ground 
state will occur. A Friedmannian expansion then takes place. The basin of attraction of 
the attractors extends to almost all of the phase space, and all initial conditions lead to the 
central Friedmann region. Since the higher-order terms do not introduce new critical points 
in the model, the qualitative picture given here should be of general validity. Notice that 
Fig. (fj) displays a model with M = 4, outside the successful range in NMC models, but 
inflationary in the NMDC model. 

5 Conclusions 

We have shown that nonminimal derivative couplings are an interesting source of new cos- 
mological dynamics. Generally speaking, their presence allows inflationary attractors where 
otherwise they would be absent. The quite narrow range of parameters that allows a viable 
inflation in NMC theories is here expanded to a semi-infinite range. Let us remark that 
the succession of different inflationary epochs may be used to break the scale-invariance of 
inflationary perturbation spectra0. In the model effectively considered, three inflationary 
phases (quasi-deSitter, power-law, quasi-deSitter) mav occur, depending on the parameter. 
This realizes a scenario of double (or triple) inflationtll without the need to introduce more 
fields. 

Before concluding, let us comment on another possible consequence of a NMDC cosmol- 
ogy. The inflationary attractors found in the previous section can be used to implement 
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a model of primordial phase transition along the scheme of extended inflation. Suppose a 
second scalar field ip with double- well self-interaction potential V^VO is present in the the- 
ory. During the inflationary stage driven by the new field may perform a first-order phase 
transition from the more energetic false-vacuum state to the less energetic true vacuum. The 
transition proceeds through quantum and thermal nucleation of true- vacuum bubbles, whose 
surface energy density contains the "latent heat" of the process. The bubbles grow convert- 
ing the background false- vacuum energy density into surface "kinetic energy" and eventually 
coalesce and percolate, unless the false-vacuum inflationary expansion is too rapid. If T is 
the nucleation rate, the (not normalized) probability that a point is still in the false- vacuum 



state at the time t isc 

p FV (t) = a\t)e-^\ (29) 



where 



lit) =4tt/3 f dt'Ta 3 (t') 
Jo 



1 du 



t' a(u) 



(30) 



where t = is the instant in which the nucleation begins. (Let us remark that ppv should be 
calculated in an Einstein frame; for what concerns the following discussion this precisation 
is of secundary importance.) If ppv — > the transition is completed. It is well known that 
in old inflation, where the expansion is exponential, the transition can never be completed; 
extended inflation cures this "graceful exit" problem by slowing down the inflation to a 
power-law. In our model the scale factor expands according to Eq. fl2"6| ) during the NMDC 
phase (the second stage) and according to Eq. (|l|) during the NMC phase (the third stage). 
It is the latter phase to be crucial. During the quasi-deSitter NMDC expansion the false 
vacuum grows so fast that the bubbles nucleated are rapidly diluted before they can coalesce, 
unless T is unrealistically large. This is the same problem one has in old inflation. However, 
during the subsequent NMC epoch the false vacuum slows down to a power-law expansion 
(it does not matter here whether inflationary or not) and the phase transition can find a 
natural exit, just as one has in extended inflation. One should also ensure that the possible 
last inflationary phase in the minimally coupled regime, which depends on the details for 
|0| — ^ of V(4>), be negligible. In this scenario, the transition occurs almost entirely during 
the NMC epoch, when the scale factor expands slowly enough. Also the bubble nucleation 
will take place mostly in the NMC epoch. Depending on the relative duration of the NMDC 
inflationary phase and the NMC phase one can have different bubble spectra. If the quasi- 
deSitter NMDC phase lasts for a large fraction of the 60 or so e-foldings one needs for the 
inflation to be successful, the bubbles produced during the subsequent power-law NMC phase 
are very small, they rapidly thermalizes (perhaps with a production of gravitational waves E|h 
and do not affect neither the large scale structure, nor the cosmic microwave background". 
If, on the contrary, the NMC phase takes, say, 55 e-foldings, then the bubbles would be 
today very large and hardly thermalized. The cosmic microwave background measurements 
and the primordial nucleosynthesis would put very stringent limits on the model, but the 
intriguing possibility that the primordial bubbles do participate in the large-scale structure 
formation arises@. 
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Figure Caption 



1. Poincare phase space for a model with m — 1,M — 3, £ = fi = — T,A = 3/2. The 
inflationary attractors are clearly seen. The trajectories start at £ = from the "north" 
and "south" poles and rapidly reach the attractors. Eventually, all trajectories fall onto 
the global stability point at = = 0, after a number of damped oscillations. The 
discrete symmetry 0, — > —0, —0, also evident in the Lagrangian, shows up. 

2. Poincare phase space for a model similar to the precedent, but with M — 4. In a NMC 
theory this model would not be inflationary. 
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